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Abstract. NLTE line formation calculations of Fe I in the solar atmosphere are extended to include weak lines in the visual 
spectrum of the Sun. Previously established atomic models are used to discriminate between different ways of treating coUi- 
sional interaction processes. As indicated by the analysis of strong Fe I lines, the influence of deviations from LTE in the solar 
atmosphere on the Fe abundance is small for all lines. To derive a common solar Fe I abundance from both strong and weak 
lines fine-tuning of the microturbulence velocity parameter and the van der Waals damping constants is required. The solar Fe I 
abundances based on all available /-values are dominated by the large scatter already found for the stronger lines. In particular 
the bulk of the data from the work of May et al. and O'Brian et al. is not adequate for accurate abundance work. Based on 
/-values measured by the Hannover and Oxford groups alone, the Fe I LTE abundances are log £fc i,© — 7.57 for the empirical 
and log£Fei,0 = 7.48 . . . 7.51 for the line-blanketed solar model. The solar Fe ionization equilibrium obtained for different 
atomic and atmospheric models rules out NLTE atomic models with a low efficiency of hydrogen collisions. At variance with 
Paper I, it is now in better agreement with laboratory Fe II /-values for all types of line-blanketed models. Our final model as- 
sumptions consistent with a single unique solar Fe abundance log epe,© ~ 7.48 . . . 7.51 calculated from NLTE line formation 
are (a) a line-blanketed solar model atmosphere, (b) an iron model atom with hydrogen collision rates 0.5 < Sh < 5 times the 
standard value to compensate for the large photoionization cross-sections, (c) a microturbulence velocity = 1-0 km s~^, (d) 
van der Waals damping parameters decreased by A log Cs = —0.10 ... — 0.15 as compared to Anstee & O'Mara's calcula- 
tions, depending on Sh, (e) Fe II /-values as published by Schnabel et al., and (f) Fe I /-values published by the Hannover and 
Oxford groups. 
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1. Introduction 

Our previous attempt to understand the formation of the iron 
spectrum in cool dwarf stars (Gehren et al. 2001 Paper I) was 
successful in isolating some of the important interaction pro- 
cesses encountered in stellar atmospheres of spectral types F 
and G. The compensating influence of (a) strong collisional 
coupling of the highly excited (> 7.3 eV) Fe I terms to the a^D 
ground state of Fe II, (b) hydrogen collision cross sections, and 
(c) photoionization from the low-excitation terms was shown 
to dominate the synthesis of line profiles and the abundances 
of solar lines. 

The lines used for the analysis were selected for strength 
because it is planned to extend the investigation to extremely 
metal-poor stars where the NLTE effects are predicted to be 
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much more important. In such stars only lines are detected 
that are strong in the Solar spectrum. The comparison of ob- 
served solar flux spectra with synthesized line profiles is thus 
hampered by all the problems usually occurring whenever line- 
broadening starts to play a role. 

The treatment of van der Waals damping had been based on 



relatively simple approximations for a long time (Unsold 1968 



Kurucz 1992), often resulting in significant underestimates of 
the damping constant. For a treatment of NLTE effects this was 
completely inacceptable, thus in Paper I we applied the quan- 
tum mechanical calculations of Anstee & O'Mara (1991 , 1995) 



without any connections. Although the results show substantial 
improvements there were still multiplets for which coiTections 
would seem adequate from profile fitting. This is not easily ex- 
plained although the calculations refer to simple LS coupling 
schemes whereas some of the upper Fe I terms involved are af- 
fected by mixing from different configurations. It appears that 
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the Anstee & O'Mara damping constants in some multiplets 
lead to line abundances that are slightly smaller than those ob- 
tained from weaker lines. 

Granular hydrodynamics are a second item that affects 
our results (Asplund et al. 2000| ). Relying on horizontally ho- 
mogeneous, plane-parallel atmospheric stratifications implies 
that dynamic movements are replaced by approximate veloc- 
ity fields, usually termed micro- and macroturbulence. For ob- 
vious reasons such an artificial replacement could depend on 
atmospheric depth as found in the empirical solar model of 
Holweger & Miiller (1974). Whereas such a stratification ^(r) 
can in principle also be constructed for other solar models, this 
is not always possible for other stars. Therefore, our fit to the 
solar Fe I line spectrum was based on a single microturbulence 
velocity ^. The values assumed for the strong lines of Paper I 
(^t = 1-00 kms^^ for the empirical and = 0.85 kms^^ for 
the line-blanketed atmospheric model) were smaller than usu- 
ally adopted for both types of model atmospheres. Thus, based 
on turbulence lines alone (lines whose equivalent widths are 
dominated by broadening due to microturbulence velocities), 
the abundances derived for both Fe II and Fe I would be slightly 
too high. 

After having examined more than 100 strong Fe I lines aris- 
ing from excitation energies between and 5 eV including 
some of the stronger turbulence lines we have found that com- 
binations of certain atomic model properties lead to acceptable 
solar flux profile fits if varying macroturbulence velocities S^t 
(Gray 1977) are applied. Due to the fact that a plane-parallel 
atmospheric model can not represent granular hydrodynamics 
with infinite accuracy, we have not tried to improve our NLTE 
profile fits beyond certain limits that are characterized by ^ 1% 
rms deviation from the observed fluxes. Yet it became clear that 
atomic models with different strengths of coUisional interaction 
led essentially to similarly good fits. This could be explained 
as a consequence of different Fe I abundances or uncertain /- 
values and van der Waals damping parameters. Unfortunately, 
the solar Fe II abundances are at least as uncertain due to sig- 
nificantly different sets of /-values. Thus the solar ionization 
equilibrium of iron could not be established because the abso- 
lute abundances were uncertain from both ends. 

As explained above part of the uncertainty remaining af- 
ter modelling the strong lines is due to line-broadening by mi- 
croturbulence and damping. Our understanding of the kinetic 
equilibrium of Fe I could therefore be considerably improved 
by extending the NLTE line formation analysis to lines that are 
substantially weaker than those of Paper I. Such lines would not 
be detected in metal-poor stars, but they would help to select 
the atomic model producing the best fit to the solar spectrum. 
Our present investigation is thus extended to a large number 
of lines with equivalent widths smaller than ^ 100mA. This 
includes lines of all degrees of excitation, although recently 
identified Rydberg transitions in the infrared with excitation 
en ergies well above 7 eV (Johansson et al. 1994, Schoenfeld et 
al. 1999) were excluded because no /-values are available. The 
following section gives a short representation of the assump- 
tions concerning both atomic and atmospheric models. Section 
3 introduces the sample of Fe I lines with results of NLTE line 
formation and profile synthesis. The last section presents our 



conclusions and a comparison with those of Paper I. We note 
in advance that the present analysis is still not able to produce 
a unique atomic model that can be applied to all kinds of stars. 
Such an investigation is left to a forthcoming paper, in which 
we will extend the analysis to a number of (mostly metal-poor) 
reference stars. 

2. Model assumptions 

2. 1. Atomic models 

Basic atomic models are the same as those of Paper I. Because 
they are described there at considerable length we will not re- 
peat the details here. The main differences between them are 
characterized by 

- the strength of the neutral hydrogen collisions, represented 
by a collision enhancement factor, S'h, which is in the 
case of no hydrogen collisions. All other cases describe 
the fact or with whi ch the collision formula propos ed by 
Drawin ( |1968| , |1969| ; see also Steenbock & Holweger [19841 ) 
is multiplied. We note that S'h ^ cxd leads to LTE. Our 
final choice resulted in S'h = 5, a value that is signifi- 
cantly greater than found previously for other atoms such 
as Al (B aumuller & Gehren |1994 |1997| ) or Mg (Zhao et 
al. 1998). Note that the role of hydrogen collisions is more 
important for Fe I than it is for Al I or Mg I, because pho- 
toionization of Fe I levels is substantially stronger than that 
of the other atoms for levels of all excitation energies; the 
large value of Sh is therefore to be considered as a com- 
pensation for the large photoionization cross-sections cal- 
culated by Bautista ( |1997[ ). 

- the treatment of the highly excited levels of Fe I. Due to 
the strong photoionization from virtually all Fe I levels the 
collisional coupling between levels above a certain limit 
^'min of excitation energy and between these levels and the 
Fe II parent terms is of critical qu ality. E lectron collisions 
are treated by the van Regemorter (1962i approximation in 
case of allowed bb collisions, by that of Allen ( 1973) for 
forbidden bb collisions, and according to Seaton's (1962) 
recipe for /j/collisions. As is obvious already from the year 
of appearance of these references collisions are the "weak 
point" of our considerations. At optical depths of the solar 
atmosphere from where most of the Fe I lines emerge, the 
resulting interaction by electron collisions is too weak to 
produce a tight coupling of the higher terms to the contin- 
uum. As a consequence, hydrogen collisions tend to result 
in a relative thermalization of only the lower Fe I terms (see 
Paper I, Figs. 6b and 6h). Thus even with strong hydrogen 
collisions (Sh > 1) only the source functions are thermal- 
ized but not the level populations or line opacities. Such a 
situation always leads to uncomfortably strong NLTE ef- 
fects in the solar spectrum. We have therefore forced ther- 
malization with respect to Fe II of all terms above £',nin, 
where different models specified E'min — 6.7, 7.0, and 7.3 
eV, respectively. In Paper I we decided to use i?min = 7.3 
eV for the final model, because that choice guaranteed that 
none of the lines and levels investigated in the solar spec- 
trum was directly affected. 
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Table 1. TH models used in the present calculations 



-0.5 

log "^5000 

Fig. 1. Photospheric solar temperature distributions of the HM 
empirical model (dashes) and the TH line-blanketed model 
(continuous curve) 



2.2. Atmospheric models 

The two plane-parallel horizontally homogeneous atmospheric 
models used in our analysis a re the semi-empirical solar model 
of Holweger & Miiller (HM, |1974| ) and our line-blanketed so- 
lar model (TH, see Paper I). Their most important difference 
with respect to line formation is the temperature stratification, 
with ThmCt) - TthCt) ~ 150 K at optical depths between 0.1 
and 1.0. The two stratifications are displayed in Fig. [l|, and the 
most important result of the temperature difference is that typ- 
ically the stronger lines are calculated with weaker line wings 
in the empirical solar model. Therefore a proper fit of Fe I line 
profiles using the HM empirical model always requires slightly 
higher damping parameters than for the TH model. 

Other important parameters of the models are those de- 
termining non-thermal spectral line core broadening. In Paper 
I we have chosen ^ = 1.00 (HM) and 0.85 kms'^ (TH), 
Srt = 2.5 (HM) and 3.2 kms^^ (TH), respectively. There is 
clear evidence that both micro- and macroturbulence vary with 
depth of line formation, however, only S^t was allowed to vary 
between '^2.0 kms~^ for some of the most saturated Doppler 
profiles and ~ 4.0 kms~^ for very weak lines. For a more re- 
aUstic analysis of both weak and strong lines in this paper we 
have added a second value of ^ = 1.00 kms^^ for the TH 
model and recalculated the non-LTE populations and line pro- 
files. No such alternative was examined for the HM model al- 
though this would probably reduce the solar Fe I abundances 
by similar amounts as for the TH model. 

The empirical HM model is used here only as a compari- 
son for abundance discussions. It had been established as a ref- 
erence for LTE conditions in the solar photosphere, and there- 
fore we have not attempted to calculate non-LTE populations 
for its temperature distribution. All the other level populations 
in this paper thus refer to the TH model for which we distin- 
guish between the (sets of) model assumptions given in Table 
|l|. Here, and Smin refer to the model atom interaction de- 
scribed in section 2.1, whereas A log Ce in the last two entries 
specifies a decrease of the damping constants with respect to 
the Anstee & O'Mara standard. This type of model leads to a 
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substantially improved fit of turbulence lines and those broad- 
ened by van der Waals damping (see below). Note that for each 
model both NLTE populations and line profiles have been re- 
calculated. When deriving solar Fe I abundances in section 3, 
some of these models are used to interpolate between different 
damping parameters (model 7 and 8). 

3. The solar weak line spectrum 

Iron is the element with probably the greatest number of lines 
visible in the solar spectrum. This is the combined result of 
a relatively high element abundance and of a very complex 
atomic configuration. In particular for Fe I nearly 1 0000 lines 
have been identified in the laboratory (Nave et al. 



1994), and 



possibly hundreds of thousands more are too weak to be de- 
tected. However, for only a small subset of these lines accurate 
/-values are known; most of them are laboratory data while 
only a subset has been derived from the solar spectrum itself. 
Our ability to identify the lines with laboratory /-values in 
the solar spectrum and calculate their solar Fe I abundances is 
therefore strongly influenced by the accuracy of the data, and it 
is this dependence that makes an analysis of the complete solar 
iron spectrum next to impossible as we will demonstrate below. 

The term "weak line" refers to all line strengths that had 
not been considered in Paper I, and it does not necessarily indi- 
cate a particularly small line strength. Thus, all lines in the list 
of Nave et al. have been examined if an /-value was available. 
Among them were only ~ 500 lines with equivalent widths be- 
low 100 mA that were not too strongly blended by other lines. 
Some of the lines retained in our sample are still blended but 
are either well-resolved or at least permit the analysis of one 
line wing. From this list we had to exclude lines in spectral re- 
gions that in the solar spectrum were overly affected by weak 
line haze and continuum uncertainties. These lie in the blue- 
green (4400 . . . 4800 A) and in the yeUow (5500 . . . 5900 A). 
The source of these spectral impurities is unknown although 
part of the blue could well be contaminated by a complicated 
pattern of Fe I autoionization transitions. Bautista's (1997) cal- 
culations show that they are there, but the accuracy of their 
wavelength positions is probably not very high. The total num- 
ber of Fe I lines including weak and strong lines was therefore 
reduced to 410, and during subsequent NLTE analyses their 
number once again shrank to the final value of 391 lines. 
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Table 2. Fe I lines in the solar flux spectrum including lines of Paper I, which have been recalculated with current model settings 
and /-values eliminating some eiTors in the previous data set. Sources of /-values and remarks are noted at the end of the table. 
NLTE models are described in the text. Equivalent widths are in mA. The complete table is only reproduced in the A&A 
version 



A [A] B[cV] logg/ 



logCg I log e(Po 1)0 (0.85) 



loge{FoI)0 (1.00) 



























LTE 


0+ 


5+ 


5- 


HM 


LTE 


5+ 


1+ 


1+ 

-0.4 


0.5+ 
-0.4 




1 


a=Do 


-z~D° 


5250 


216 





121 


-4 


94 


e 


-32 


051 


7.57 


7.71 


7.61 


7.74 


7.66 


7.49 


7.53 


7.55 


7.56 


7.59 


acf 


71.2 


1 


a^Di 




5225 


533 





110 


-4 


79 


e 


-32 


052 


7.60 


7.73 


7.64 


7.77 


7.69 


7.51 


7.55 


7.57 


7.59 


7.60 


af 


76.1 


1 


a=D2 


-z^U° 


5247 


057 





087 


-4 


95 


e 


-32 


057 


7.56 


7.72 


7.59 


7.75 


7.65 


748 


7.53 


7.55 


7.57 


7.58 


af 


68.7 


1 


a^D4 


-z^D° 


5166 


282 





000 


-4 


20 


e 


-32 


070 


7.51 


7.66 


7.64 


7.71 


7.69 


7.42 


7.51 


7.51 


7.52 


7.53 


cdfh 


109.0 


2 






4445 


480 





087 


-5 


44 


e 


-31 


996 


7.54 


7.65 


7.54 


7.66 


7.65 


7.50 


7.52 


7.53 


7.54 


7.55 




43.3 


2 


a=D3 


-z^F° 


4427 


309 





052 


-2 


92 


a 


-32 


000 


7.58 


7.66 


7.58 


7.65 


7.73 


7.54 


7.56 


7.56 


7.69 


7.70 


dci 


199.1 


2 


a^D4 




4347 


237 





000 


-5 


50 


e 


-32 


006 


7.54 


7.65 


7.55 


7.66 


7.66 


7.49 


7.51 


7.53 


7.53 


7.54 


ac 


43.6 


3 


a^Di 


-z^P2 


4232 


720 





110 


-4 


93 


e 


-31 


968 


7.47 


7.57 


7.49 


7.59 


7.59 


7.40 


7.44 


7.45 


7.46 


7.47 




60.2 


13 


a^Fl 


-z^F? 


6625 


026 


1 


Oil 


— 5 


35 


m 


-31 


927 


7.52 


7.64 


7.53 


7.65 


7.63 


7.51 


7.53 


7.54 


7.54 


7.55 


iig 


16.4 


13 


a^Fa 


-z^F5 


6574 


233 





990 


— 5 


02 


a 


-32 


400 


7.56 


7.69 


7.57 


7.69 


7.68 


7.54 


7.56 


7.58 


7.58 


7.60 




29.7 


13 


a=F3 


-z^F° 


6498 


945 





958 


-4 


70 


f 


-31 


934 


7.58 


7.70 


7.59 


7.73 


7.67 


7.54 


7.56 


7.58 


7.59 


7.61 


ag 


48.1 


13 


a^F4 




6400 


323 





915 


-4 


32 


a 


-31 


939 


7.49 


7.60 


7.51 


7.63 


7.56 


7.44 


7.45 


748 


7.49 


7.51 


bf 


65.5 


13 


a=F5 


-z^F? 


6280 


620 





859 


-4 


39 


f 


-31 


947 


7.55 


7.71 


7.60 


7.75 


7.66 


7.52 


7.52 


7.52 


7.53 


7.53 


dc 


68.4 


14 


a^F4 




6120 


250 





915 


-5 


95 


m 


-31 


926 


7.53 


7.64 


7.54 


7.65 


7.69 


7.53 


7.53 


7.53 


7.53 


7.53 


gc 


5.6 


15 


a=F2 


-z^D° 


5405 


775 





990 


-1 


88 


P 


-31 


870 


7.50 


7.52 


7.51 


7.60 


7.64 


7.47 


7.49 


7.50 


7.64 


7.64 


e 


271.5 


15 


a=F3 


-z^D" 


5371 


489 





958 


-1 


65 


n 


-31 


870 


741 


7.44 


7.43 


7.52 


7.57 


7.40 


7.40 


7.40 


7.53 


7.53 




306.3 


15 


a^F4 


-z^D° 


5328 


038 





915 


-1 


47 


n 


-31 


880 


7.46 


7.52 


7.47 


7.56 


7.62 


7.44 


7.46 


7.47 


7.61 


7.61 


e 


397.9 


15 


a=F4 


-z^I)° 


5397 


128 





915 


-1 


99 


n 


-31 


880 


7.47 


7.51 


7.48 


7.56 


7.64 


7.46 


7.46 


7.47 


7.63 


7.63 


i 


241.7 


15 


a=F5 




5269 


537 





859 


-1 


32 


n 


-31 


890 


7.44 


7.51 


7.44 


7.54 


7.61 


7.43 


7.45 


7.45 


7.61 


7.61 




501.5 


34 


a^Fa 


-z^Fj 


6851 


640 


1 


608 


-5 


32 


b 


-31 


786 


7.46 


7.58 


7.47 


7.58 


7.60 


7.46 


7.47 


7.49 


7.49 


7.50 


acf 


3.9 


34 


a^Fa 


-z^F" 


6739 


540 


1 


557 


-4 


79 


P 


-31 


795 


7.37 


7.49 


7.38 


7.52 


7.48 


7.37 


7.38 


7.40 


7.40 


7.41 


acf 


12.3 


34 


a^Fi 


-z5f° 


6581 


220 


1 


485 


-4 


68 


P 


-31 


806 


7.42 


7.54 


7.43 


7.55 


7.54 


741 


7.42 


7.44 


7.44 


7.45 


ac 


20.5 


34 


a^Fi 


-z^F° 


6710 


310 


1 


485 


-4 


88 


m 


-31 


813 


7.54 


7.65 


7.55 


7.66 


7.64 


7.53 


7.54 


7.55 


7.55 


7.56 


acf 


16.2 


36 


a^Fs 


-z^F° 


5216 


274 


1 


608 


-2 


15 





-31 


670 


7.45 


7.45 


7.48 


7.59 


7.53 


7.35 


7.39 


741 


7.51 


7.51 


dcfh 


130.3 


36 


a^Fa 


-z3f° 


5194 


941 


1 


557 


-2 


09 





-31 


680 


7.45 


7.46 


7.45 


7.55 


7.55 


7.36 


7.38 


7.39 


7.50 


7.50 


ch 


129.6 


38 


a^F2 




4798 


734 


1 


608 


-4 


25 


b 


-31 


612 


7.59 


7.70 


7.60 


7.71 


7.71 


7.58 


7.59 


7.60 


7.60 


7.61 


gh 


34.3 


38 


a^Fa 




4772 


820 


1 


557 


-2 


90 


a 


-31 


627 


7.68 


7.81 


7.72 


7.86 


7.81 


7.61 


7.65 


7.66 


7.72 


7.73 


dg 


93.7 


41 


a^F3 


-z5g° 


4404 


750 


1 


557 


-0 


10 


P 


-31 


560 


7.40 


7.48 


741 


7.50 


7.57 


7.39 


7.41 


741 


7.54 


7.55 


h 


786.0 


41 


a^Fi 




4383 


545 


1 


485 





20 





-31 


580 


7.39 


7.49 


7.42 


7.48 


7.58 


7.39 


7.41 


741 


7.55 


7.56 


h 


1345.6 


42 


a^F4 


-z3g° 


4147 


669 


1 


485 


-2 


10 





-31 


520 


7.47 


7.51 


7.48 


7.58 


7.58 


741 


7.41 


741 


7.51 


7.51 


ch 


131.2 


42 


a^'Fi 


-z3g§ 


4271 


760 


1 


485 


-0 


16 





-31 


550 


7.37 


7.44 


7.38 


7.46 


7.55 


7.36 


7.36 


7.36 


7.51 


7.51 


h 


846.7 


43 


a^Fa 




4071 


738 


1 


608 


-0 


02 





-31 


440 


7.33 


7.41 


7.33 


7.41 


7.51 


7.31 


7.31 


7.31 


7.47 


7.47 


h 


860.9 


43 


a^Fs 


-!/^F° 


4063 


594 


1 


557 





06 


a 


-31 


470 


7.31 


7.45 


7.35 


7.46 


7.55 


7.35 


7.35 


7.36 


7.50 


7.50 


hj 


900.8 


43 


a^'Fi 


-!/='F2 


4045 


812 


1 


485 





28 





-31 


490 


7.36 


7.43 


7.38 


7.46 


7.55 


7.36 


7.36 


7..36 


7.51 


7.51 


hj 


1250.5 


62 






6297 


800 


2 


223 


-2 


73 


q 


-31 


565 


7.57 


7.67 


7.58 


7.71 


7.64 


748 


7.50 


7.52 


7.53 


7.55 


h 


75.3 


62 


a^P3 




6151 


620 


2 


176 


-3 


27 


q 


-31 


569 


7.51 


7.62 


7.51 


7.61 


7.55 


7.45 


7.46 


7.47 


7.48 


7.50 


h 


51.3 


63 


a^Pl 




6015 


250 


2 


223 


-4 


68 


m 


-31 


539 


7.52 


7.62 


7.53 


7.63 


7.65 


7.52 


7.53 


7.55 


7.54 


7.55 




4.5 


64 


a^Pl 


-z3p° 


6082 


720 


2 


223 


-3 


59 


P 


-31 


545 


7.53 


7.63 


7.54 


7.65 


7.61 


7.51 


7.52 


7.53 


7.54 


7.55 


r 


35.8 


64 


a^Pl 


-z3po 


6240 


660 


2 


223 


-3 


23 


P 


-31 


560 


7.48 


7.57 


7.49 


7.61 


7.55 


7.43 


7.44 


7.45 


7.47 


7.48 




50.3 


66 


a^Pl 




5198 


711 


2 


223 


-2 


14 


g 


-31 


440 


7.58 


7.65 


7.60 


7.70 


7.63 


7.51 


7.53 


7.54 


7.56 


7.57 


dfk 


103.9 


66 


a = P2 


-!/^P? 


5079 


223 


2 


198 


-2 


07 


g 


-31 


430 


7.55 


7.64 


7.59 


7.68 


7.67 


7.50 


7.52 


7.53 


7.57 


7.59 


bfk 


107.5 


66 


a^Pa 


-!/'P5 


5145 


099 


2 


198 


-2 


88 


a 


-31 


439 


7.26 


7.36 


7.27 


7.38 


7.34 


7.22 


7.24 


7.25 


7.27 


7.28 


d 


54.3 


66 


a^Pa 


-!/'Pl 


5250 


646 


2 


198 


-2 


18 


a 


-31 


460 


7.69 


7.75 


7.72 


7.82 


7.76 


7.59 


7.62 


7.65 


7.72 


7.73 


ch 


108. 1 


68 


a^Pl 




4447 


717 


2 


176 


-1 


34 


g 


-31 


270 


7.61 


7.70 


7.64 


7.73 


7.79 


7.58 


7.61 


7.63 


7.74 


7.76 


d 


185.4 


68 


a^Pa 




4494 


563 


2 


198 


-1 


14 


g 


-31 


300 


7.48 


7.53 


7.50 


7.59 


7.61 


7.45 


7.47 


748 


7.60 


7.60 


df 


206.8 


69 


a'Spa 




4447 


130 


2 


198 


-2 


73 


a 


-31 


280 


7.63 


7.72 


7.64 


7.74 


7.68 


7.55 


7.57 


7.59 


7.61 


7.63 


d 


66.4 


69 


a^P3 


-s/'^PS 


4442 


840 


2 


176 


-2 


79 


g 


-31 


290 


7.55 


7.64 


7.55 


7.64 


7.64 


7.51 


7.51 


7.51 


7.53 


7.54 


bl 


64.5 


71 


a^P3 




4282 


402 


2 


176 


— 


78 


a 


— 31 


240 


7.17 


7.23 


7.18 


7.26 


7.33 


7.15 


7.15 


7.15 


7.28 


7.29 


bl 


193.8 


109 


a^Pa 




6392 


543 


2 


279 


-4 


03 


m 


-31 


553 


7.57 


7.68 


7.58 


7.69 


7.66 


7.56 


7.57 


7.58 


7.58 


7.59 


g 


19.0 


109 


a^Pa 




6608 


030 


2 


279 


-4 


03 


b 


-31 


570 


7.56 


7.66 


7.56 


7.67 


7.65 


7.55 


7.56 


7.57 


7.57 


7.58 


g 


18.0 


111 


a^Po 


n 


6978 


850 


2 


484 


-2 


48 


P 


-31 


523 


7.59 


7.64 


7.63 


7.73 


7.62 


7.53 


7.54 


7.57 


7.61 


7.62 


af 


78.1 


111 


a^Pl 


-z3p° 


6663 


450 


2 


424 


-2 


45 


P 


-31 


521 


7.54 


7.60 


7.58 


7.70 


7.58 


748 


7.50 


7.52 


7.56 


7.57 


afh 


80.8 


111 


a^Pl 




6750 


150 


2 


424 


-2 


61 


P 


-31 


528 


7.60 


7.67 


7.63 


7.73 


7.63 


7.53 


7.54 


7.56 


7.60 


7.61 


af 


76.7 


111 


a^Pa 


- z-^P^ 


6421 


350 


2 


279 


-1 


95 


P 


-31 


560 


7.45 


7.48 


7.49 


7.60 


7.52 


7.38 


7.43 


7.44 


7.52 


7.52 


bth 


110.0 


113 


a^Pl 




5678 


600 


2 


424 


-4 


67 


m 


-31 


426 


7.51 


7.61 


7.51 


7.61 


7.65 


7.51 


7.52 


7.53 


7.53 


7.54 


cd 


3.2 


113 


a^Pa 




5436 


590 


2 


279 


-2 


96 


a 


-31 


451 


7.19 


7.29 


7.20 


7.31 


7.27 


7.16 


7.17 


7.18 


7.19 


7.20 


b 


45.6 


114 


a^Pl 




5141 


739 


2 


424 


-1 


96 


P 


-31 


350 


7.37 


7.42 


7.38 


7.49 


7.41 


7.36 


7.29 


7.30 


7.36 


7.37 


dh 


88.1 


114 


a^Pa 




4924 


769 


2 


279 


-2 


24 


q 


-31 


370 


7.71 


7.74 


7.73 


7.83 


7.78 


7.64 


7.64 


7.65 


7.70 


7.72 


df 


97.6 


114 


a^Pa 




5049 


819 


2 


279 


-1 


33 


q 


-31 


390 


7.51 


7.53 


7.52 


7.62 


7.65 


7.46 


7.48 


7.49 


7.60 


7.62 


df 


164.9 


115 


a^Pa 


Da 


4574 


720 


2 


279 


-2 


97 


b 


-31 


278 


7.65 


7.77 


7.68 


7.78 


7.73 


7.60 


7.61 


7.62 


7.65 


7.67 


a 


59.8 


116 


a^Pa 


-z^S° 


4439 


880 


2 


279 


-3 


00 


g 


-31 


237 


7.52 


7.62 


7.53 


7.63 


7.61 


7.48 


7.49 


7.50 


7.51 


7.52 


a 


52.9 


152 


z'''D° 


-e''Da 


4233 


602 


2 


482 


-0 


60 


g 


-30 


640 


741 


7.46 


7.39 


7.46 


7.55 


7.37 


7.37 


7.37 


7.51 


7.51 


cf 


278.4 


152 


z^-D° 


-e''D3 


4250 


119 


2 


469 


-0 


41 


g 


-30 


660 


7.45 


7.53 


7.46 


7.54 


7.63 


7.44 


7.46 


7.46 


7.61 


7.62 


ade 


355.4 


152 




-e''Da 


4187 


039 


2 


449 


-0 


55 


g 


-30 


640 


7.39 


7.47 


7.40 


7.49 


7.55 


7.38 


7.39 


741 


7.54 


7.54 


ade 


297.2 


152 


z^I)° 


-e''D3 


4222 


213 


2 


449 


-0 


97 


g 


-30 


650 


741 


7.51 


7.44 


7.53 


7.56 


741 


7.42 


7.42 


7.56 


7.56 


ade 


198.5 



(1976), (g) Blackwell et al. (|1982 a 
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Line synthesis remarks: (a) no blend, no asymmetry, (b) resolved blend(s), (c) continuum adjusted, (d) unresolved blend(s), (e) core 
asymmetry, (f) blue and red wing deficit, (g) only red wing deficit, (h) core too wide, (i) core too narrow, (j) core too deep, (k) core too shallow, 
(1) all faint lines included 
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One of the more surprising results of this evaluation of the 
solar Fe I spectrum is that the number of truly weak lines with 
both an acceptable spectral environment and laboratory /-value 
is so small. This is the case for lines in a range of solar equiva- 
lent widths from 3 to 30 mA. This has also been noticed among 
others by Rutten & van der Zalm (1984). If laboratory analy- 
ses were extended into the near infrared the line list could be 
greatly extended because of decreasing blend problems. The 
blue and near-ultraviolet spectral regions have been ignored 
here because of the problems localizing the continuum below 
4200 A. 

3. 1. Spectrum synthesis 

The final set of Unes is reproduced in Table ^ together with all 
relevant data. The sources of the /-values as well as the re- 
marks in the second last column are noted at the end of the ta- 
ble. The damp ing constan ts are calculated according to Anstee 
& O'Mara's ( 1991 , 1995 ) theory as in Paper I, and they are 
given here in terms of van der Waals damping constants. The 
equivalent widths in the last column are integrated on the basis 
of the best synthetic fit of the solar ^mjc profile. We emphasize 
that they are not used for the line analysis which is solely based 
on profile fits. Rather, they are derived from the theoretical pro- 
file after the final profile fitting procedure. Their accuracy is 
low, which is uncritical since they are used for graphical pur- 
poses only. 

3.1 .1 . Oscillator strengths 

In order to determine abundance ratios in spectral lines of stars 
other than the Sun it is often sufficient to know the product 
gfsQ, which can be obtained in the solar flux spectrum with 
no particular knowledge of the /-value. Were it not for con- 
sistency and identification checks and for the determination of 
the solar iron abundance itself, no oscillator strengths would 
be needed. Such consistency checks include the specification 
of broadening parameters such as microturbulence and damp- 
ing constants, because both can to a certain degree replace 
abundances or oscillator strengths. Therefore a critical analy- 
sis of the /-values is necessary. As mentioned above, oscillator 
strengths available for Fe I lines come from essentially three 
different methods: 

- Theory has made important progress in the last 20 years. 
This is not only seen in the bf cross-sections we used in our 
kinetic equilibrium calculations but also in a virtually com- 
plete set of calculated /-values made available by Kurucz 
(1992 1. The main obstacle in using these data lies in the 
problem of estimating their accuracy. Therefore we have 
used Kurucz' data for their original statistical purpose com- 
puting particle interaction via the statistical equilibrium 
equations whenever laboratory /-values were missing. But 
we have not applied his /-values during subsequent spec- 
trum synthesis. 

- Laboratory methods have made some progress, too, and the 
number of laboratory /-values is steadily increasing. It is 
this origin we have put most confidence in, although the 



reliability of the various sources as judged from their ability 
to fit the solar flux spectrum is surprisingly different as we 
will show below. 

The inversion method, i.e. measuring solar /-values by 
synthesis of solar equivalent widths, has become a popu- 
lar method to fill the missing data gap in the Fe I line list. 
Whether based on equivalent widths or line profiles, this 
method always reproduces an assumed abundance scale. 
This is - in most cases - the meteoritic Fel abundance, 
sometimes it is tied to some otherwise established solar iron 
abundanc e such as that of the Oxford group (cf. Gurtovenko 
& Kostik 1981). It never carries information about the os- 



cillator strength itself. 

Of these three methods we have applied only the results of the 
last two methods to spectrum synthesis, and in solar abundance 
determinations we confine our sample to those lines for which 
laboratory /-values are available. 

3.1.2. Line broadening 

During the analysis of the strong Fe I lines we have discussed 
collision broadening at some length in Paper I, where it was 
documented that the results of the broadening theory of Anstee 
& O'Mara ( |l99l|Jl995| ) provided the necessary adjustment be- 
tween weak and strong line abundances at least in a qualitative 
way. We have followed this approach in the present investiga- 
tion, replaced the old van der Waals damping constants by the 
new collision parameters, however, staying aware of the depen- 
dence of abundance analyses upon atmospheric models. Thus, 
the empirical model of Holweger & Miiller (1974) requires sig- 
nificantly higher damping constants than our theoretical model 
in order to fit solar strong line profiles with the same abundance 
as the weak lines. In fact, we have added to our sample of NLTE 
models two more items with reduced damping constants in or- 
der to explore their influence on the mean Fe I abundance. We 
come back to this point in subsection 3.2.2. 

The introduction of weak lines, among them many lines 
broadened by microturbulence, has considerably enhanced our 
possibility to judge the solar line spectrum and the necessary 
atomic data. So the present analysis required an extension of 
the parameter space covered by non-thermal motions to put 
both weak and turbulence lines on a common abundance level. 
In fact, iiTespective of the source of /-values, lines between 
50 and 120 mA tend to require systematically higher abun- 
dances than weak or very strong lines if the value of Paper 
I, ^ = 0.85 kms^^ was used. We introduced a second mean 
value of ^ = 1.00 kms~^ which seems more appropriate for 
our present investigation. Note that this value has only limited 
influence on the strong lines, so our former results stay essen- 
tially unchanged. 

the details of turbu- 



As will be shown in subsection 3.1.3 



lent line broadening are still unsatisfactory for a number of 
medium-strong lines. Whereas all weak lines with equivalent 
widths below W\ ^ 70 mA and most of the very strong lines 
are well represented by the synthetic line profiles, some lines 
around Wx ^ 70 . . . 120 mA are not reproduced by any choice 
of model parameters. This was noticed already in Paper I when 
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Fig. 2. Abundance differences between lines synthesized in our 
plane-parallel LTE (TH) model and those obtained from a hy- 
drodynamical solar model of Asplund et al. ( 2000| ). Lines that 
were synthesized in our plane-parallel model with continuum 
adjustment are drawn as open circles 



trying to fit Fe II multiplet 42 or Fe I multiplets 1 or 36. The 
present selection of Fe I lines includes quite a lot of such lines 
that seem to document the ultimate difference between plane- 
parallel and hydrodynamical models. Following this difference 
it is interesting to compare the results of the two completely 
different model realizations of non-thermal motions. Therefore 



the results of Asplund et al. ( 200C ) have been confronted with 
our data in Fig. ^ 

It is true that the mean abundance of the 49 lines in com- 
mon is different by A log e ~ 0.05 (or even slightly more for 
turbulence lines), and this could be interpreted as the differ- 
ence between plane-parallel and hydrodynamical models. But 
a closer view reveals that most of the weaker lines belong to 
a category that requires some continuum adj ustme nt with re- 
spect to the solar flux atlas of Kurucz et al. ( 1984 ). There are 
some spectral regions that suffer from unknown continuum de- 
pressions, and whenever such an adjustment was used in our 
calculations, the abundance differences between our respective 
models shrank to a mean A lege ~ 0.03, more probably near 
the true difference between the models. It is interesting in this 
respect that the bulk of turbulence line abundances between 60 
and 90 mA is systematically higher than those calculated from 
the hydrodynamical model. This is also found in our own data 
when strong lines and turbulence lines are compared, and it 
would mean that exactly this type of lines is not particularly 
well synthesized by plane-parallel models. 

We emphasize, however, that a single value for the micro- 
turbulence velocity cannot be assumed to reproduce all types of 
core saturation found in turbulence lines. Our simple approxi- 
mation is inconsistent in that it ignores the corresponding varia- 
tions found and accepted for the macroturbulence velocity, and 
a free fit of the ^ parameter for each line profile would have pro- 
duc ed sli ghtly improved results. Comparison with Asplund et 
al. (2000) finally shows that both weak and strong lines are not 
strongly affected by dynamic processes, which means that the 
conventional replacement of laminar flow patterns by a micro- 
/macroturbulence approach is still surprisingly valid. 



3.1 .3. Line profiles and equivalent widths 

The overwhelming majority of publications is devoted to the in- 
vestigation of equivalent widths which is mostly due to the easy 
access to such data in the literature. The critical examination of 
line profiles instead makes available an increased amount of 
information about line formation and stellar atmospheric con- 
ditions. Our present work on NLTE effects in Fe I lines is based 
on roughly 4000 line profiles, and their evaluation is coded in 
a very coarse set of remarks in Table ^. Such remarks combine 
the average profile properties of all models for a particular line, 
and the following description will show only typical properties. 

Very weak lines (W\ < 10 mA): 

Only 10% of the total sample consist of very weak lines. Most 
of them could be selected to be free from known blends, but 
only 10 of them were unaffected by problems with continuum 
adjustment. It is this latter quality that makes the analysis of 
very weak lines so ambiguous. This can be seen in Fig. ^ where 
the LTE profile fits for two lines are shown. Continuum adjust- 
ment is by far not always as small as 0.5% as it is for the line in 
Mult 1 109, and ignoring it may lead to abundances higher by 
up to 0.15 dex in single cases. 

It is no straightforward procedure to decide which lines to 
submit to continuum adjustment, because this requires a look at 
the whole spectral region. Consequently, we have adjusted the 
atlas continuum only if there is a continuum depression over 
at least 10 A. In some cases we tried to synthesize faint back- 
ground lines in order to estimate their influence on the con- 
tinuum position. While weak lines should be least affected by 
broadening and therefore yield most reliable abundances, the 
continuum placement destroys a substantial part of this argu- 
mentation. 

Weak lines (10 < M^a < 60 mA): 

These lines constitute the majority of the sample with more 
than half in this range of equivalent widths. Up to 30 mA the 
lines do not depend significantly upon microturbulence, but 
their abundance change increases to —0.03 per 0.1 kms~^ at 
60 mA. A number of weak lines that are fairly representative 
of our sample is reproduced in Fig. ^ together with LTE profile 
fits for both the HM and TH models. They are shown in partic- 
ular to demonstrate the abundance differences between the two 
models. It should be mentioned here that this subsample of Fe I 
lines produces by far the best profile fits, followed by the strong 
lines, the very weak lines, and the turbulence lines, in order of 
decreasing fit quality. The profiles of the weak lines are not 
dictated by core saturation or line wing broadening but, nearly 
exclusively, by external line broadening due to solar rotation 
and macroturbulence. As is the case for some of the very weak 
lines, some weaker lines in Fig. ^ require a high macroturbu- 
lence of S > 4 km s^^ in order to adjust the wings. 

We note that the quality of the profile fit is the same for both 
atmospheric models, irrespective of the abundance differences. 
Thus most of the very weak and weak lines show a systematic 
abundance difference of A(log£Fc,o)HM-TH = 0.06 . . . 0.12 
(see below). As with the very weak lines, there is also no prob- 
lem when fitting the profiles of the weak lines with different 
NLTE models (not shown in Fig. 0). However, the kinematic 
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Fig. 4. Profiles of weak lines (10 < Wx < 60 mA) of Fe i in the solar flux spectrum (filled circles). Synthetic profile fits are 
LTE and HM ( ) or TH ( ) atmospheres. Fit parameters are indicated 
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Fig. 6. Profiles of turbulence lines (60 < W\ < 110 mA) of Fe I in the solar flux spectrum (filled circles). Synthetic profile fits 
are asain for LTE and HM { ^ or TH { ^ atmosoheres. Fit oarameters are indicated 
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0.96 



LTE (TH) fit with ^ = 1 .0 km/s, 2 = 3.1 km/s 




Fel 1142 (9.6 mA) 
log e(Fe) = 7.58 



6033.9 



- LTE (TH) fit with ^ = 1.0 km/s, E = 3.6 km/s [continuum adjusted, black dots] 



1.00 




Fel 1109 (8.6mA) 
log s(Fe) = 7.58 



Fig. 3. Typical problems with the adjustment of the local so- 
lar flux atlas continuum. Top: Weak line in Mult 1 142 with no 
changes of the local continuum necessary. Bottom: Mult 1109 
shown with and without continuum adjustment. The original 
atlas spectrum is reproduced with open circles and fitted by the 
grey curve with logepc = 7.63 and S = 4.0 kms^^ 



properties of all lines with equivalent widths below 100 mA are 
reproduced in a number of profiles that show systematic bisec- 
tor curvature and a red line wing deficit. An even more critical 
inspection of some of the profiles reveals synthetic line cores 
fliat tend to be too broad even for ^ = 0.85 (TH) or 1.00 (HM) 
kms~^, respectively. This is evident in particular for lines that 
are formed further up in the atmosphere, and - together with 
the red wing asymmetries - it clearly documents the pitfalls of 
static atmospheric models. Some of the weak lines are also af- 
fected by a bad definition of the local continuum, which either 
lead to a removal of a significant number of lines originally se- 
lected or ended in a multi-line synthesis with a number of faint 
background lines included. Such results are not given too much 
weight in the abundance analysis. 

Turbulence lines (60 < M^a < HO mA): 

Roughly 20% of our sample are strong enough for core satu- 
ration and are therefore shaped by the value of the microtur- 
bulence parameter. Naturally, a static model atmosphere repro- 
duces such lines only in an approximative way. This is seen 
in Fig. H where a number of such lines and their synthetic fits 
are presented. Most of these fits require substantially smaller 



Fel 66 



0.85 / LTE: TH / |Fe/H| = 0.18 / S = 3.0 
1 .00 / LTE: HM / IFe/H] = 0.25 / S = 2.4 



5250.50 



5250.60 



5250.70 



5250.80 



Fig. 5. LTE profiles of Fel 66, 5250.646 A. Models are as in 
Figs. I and|. Additionally, a TH LTE model with ^ = 1.0 
kms^^ is plotted for comparison (dashes). The deep profiles 
are uncorrected for rotation and macroturbulence, the "v"- 
shaped profiles include external broadening 

values of the macroturbulence velocity S, but even then the 
synthetic core profiles are often too broad and too shallow. 

In contrast to weaker lines for which the fit with synthetic 
profiles can be made nearly as accurate as desired, the fit of 
turbulence lines with a plane-parallel atmospheric model has 
its natural limitations which are explained by the velocity dif- 
ferences necessary to fit the innermost core and the wings si- 
multaneously. Thus, in principle the saturated core seems to 
require relatively small velocity fields, whereas the opposite is 
required for the wings, a modulation that roughly represents 
the hydrodynamic equation of continuity. The microturbulence 
values used in the LTE models of Fig. ^ have in fact been cho- 
sen so as to fit the line core width. Using even larger values 
as would be indicated by comparison with weak and strong 
lines does not improve the profile fits although it may help 
to minimize the overall abundance scatter. Fig. || emphasizes 
the difference in core saturation between the two model atmo- 
sphere types (HM and TH). Due to the temperature differences 
between the atmospheric models profiles synthesized from the 
HM model always require a smaller macroturbulence to fit the 
very line core than do the LTE or NLTE profiles based on the 
TH model. 

We note that turbulence velocity gradients introduced 
within the scope of static plane-parallel models do not improve 
the profile fits either. The kinematic fine-tuning of the turbu- 
lence lines thus will stay the exclusive domain of granular hy- 
drodynamics. 

Again, as with the weaker lines, LTE and NLTE models 
both tend to produce similar profile fits for the turbulence lines 
provided that the abundances are correspondingly adjusted. 
This is a direct consequence of the source function thermaliza- 
tion inherent to our NLTE modelling. As can be seen in Table ||, 
lines with equivalent widths around 100 mA display an abun- 
dance spread of ~ 0.2 dex among different LTE and/or NLTE 
models. 

The profiles of the stronger Fe I lines (W\ > 1 10 mA) have 
been discussed in Paper I. It is therefore sufficient to repeat 
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Fig. 7. LTE profiles of Fe I 1 197, 6726.670 A, computed with 
the HM model atmosphere displaying the sensitivity of turbu- 
lence lines with respect to abundance changes 

here, that simultaneous fits of line cores and damping wings are 
only obtained outside the range of the inner wings (±0.1 . . . 0.4 

A). 

3.2. Abundances 

Our investigation of NLTE excitation and ionization in the so- 
lar photosphere would not be complete without mentioning the 
solar Fe I abundance problem. Since there exists quite a number 
of publications on the "true " solar Fe I abu ndance (e.g. Bi emont 
et al. |l99li B lackw ell et al. |l9954|l995b[ Hol weger et al. |1995 
Kostik et al 



1996 



Grevesse & Sauval 1999 ), we will not en- 
ter into details but simply give our judgement according to the 
large number of lines of all strengths examined with reference 
to complete profile information (but ignoring their center-to- 
limb variation) and an exhaustive range of NLTE models. 

Current analyses tend to put their results into perspective by 
denoting the differences between photospheric and meteoritic 
Fe I abundances. The latter has been known for many years now 
(Anders & Grevesse 1989), logepci,© — 7.51. Photospheric 



abundance determinations, however, range from logeFci.o — 
7.42 (Schnabel et al. |l999l Fe ll) to 7.67 (Blackwell et al. 



1995a , Fel). As was pointed out by Kostik et al. (1996) and 
later iterated by Grevesse & Sauval ( 1999| ), the discrepancy be- 



tween different groups of researchers depends on a number of 
different methods and data sets the influences of which are not 
always easily disentangled. 

A few problems have already been discussed above, in par- 
ticular the important influence of selecting a local spectral con- 
tinuum. Other problems arise when determining abundances 
based on measurements of equivalent widths. Thus, Meylan et 
al. (1993) have used Voigt profile fits to reproduce their ob- 
served Fe I lines. Their results differ systematically from those 
of other methods produced either by planimeter measurements 
or - as in our case - from full line profile synthesis. This is an 
important source of systematic eiTors because anything but fit- 
ting synthesized profiles requires an estimate of the line wing 
extension that is often - and always systematically - neglect- 
ing a weak line haze. One of the more moderate examples is 



reproduced in Fig. ^. For this line Meylan et al. (1993) list an 
equivalent width of 53.6 mA, obtained from their Voigt profile 
fit. Our synthesis reproduces the observed solar flux spectrum 
with no continuum adjustment applying an Fe I abundance of 
log£Fci,0 = 7.55, whereas their equivalent width requires an 
abundance ^0.1 dex higher than ours. More importantly, their 
equivalent width does not fit the observed profile. Other turbu- 
lence lines listed by Meylan et al. show even larger discrepan- 
cies up to 0.3 dex! Therefore it is not surprising that - using 
the /-values published in that paper - we derive a mean so- 
lar abundance of logepci,© = 7.25. Altogether, at this stage of 
analyzing the solar Fe I abundance we ignore solar /-values be- 
cause they would not add to abundance information, since their 
determination requires the input of a mean abundance value. 

3.2.1 . Sources of oscillator strengths 



Except fo r the r esults of Meylan et al. ( 1993 ) and Gurtovenko 
& Kostik ( 198 1 ) Table contains only references to laboratory 



/-values that cover more than 80% of the lines. Among them 
we find essentially four different sets of data, 

- The laser-induced fluorescence measurements of O' Brian 



et al. ( |1991| ), 

- /-values obtained from stabilized arc-emission by May et 
al. ( |1974| ), 

- Observations of stabilized furnace absorption by the 



Oxford group of Blackwell et al. ( 1976 1979^ , 1979b 



198C 



1982E, 1982b) 



Hollow-cathode and laser-induced fluorescence measure- 
ments performed by the Hannover group of Bard et al. 
(|1991|,|1994|) 



The rest of the sources is not very important for our investiga- 
tion. The results listed in Table || refer to a broad selection of 
methods which have been repeatedly discussed (see Holweger 
et al. | 1995[ Kostik et al. |1996| or Grevesse & Sauval |1999[ ). 



We start with a plain characterization of the abundance results 
obtained with the different sets of /-values. 

The top frame of Fig. [s] shows LTE abundance results ob- 
tained from the HM e mpiri cal model atmos phere using the 
data of O'Brian et al. ( |l99l| ) and May et al. ( |l974 , whereas 
the bottom frame of Fig. |8j displays the results for the oscilla- 
tor strengths determined by the Oxford and Hannover groups. 
While the proper choice of models and parameters is discussed 
in the following subsection, it is already evident here that the 
two frames harbour sources of different quality. Thus, the /- 
values of O'Brian et al. or May et al. lead to approximately 
twice the r.m.s. scatter of the solar abundances as compared 
with the results derived from the /-values of the Oxford and 
Hannover groups. The May et al. abundances are also system- 
atically higher than the mean. 

The /-values of O'Brian et al. and those of Bard and Kock 



( |1994| ) are on the same absolute scale since both have used 
very similar measurements and normalization procedures. In 
fact. Fig. 3 in Bard & Kock shows a negligible difference of 
the coiTesponding /-values for the lines in common, although 
the strong scatter is confirmed. What makes the O'Brian et al. 
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sample so suspicious is the occurrence of abundance differ- 
ences between lines in a common multiplet. An extreme case 
is Mult 66, where our results for A5145.099 and A5250.646 
lead to logepci — 7.34 and 7.76, respectively. There are also 
other lines such as A4798.267 and A4735.845 of Mult 1042 
with log£Foi = 7.35 and 7.81, respectively. 

There is no simple explanation why the oscilla tor strengths 
of May et al. and those of Bard and Kock (1994) lead to dif- 
ferent abundances. The data used in our analysis are those in 
Fuhr et al. (1988), which had been renormalized to the scale 



of the Oxford measurements. Most of the corrected May et al. 
/-values are therefore 0.1 dex smaller than the original data. 
Based on the original paper, the May et al. abundances thus 
would be 0. 1 dex smaller. While this accounts for half of the 
difference between the two groups, there remains another 0. 1 
dex difference which is not seen in Fig. 2 of Bard & Kock. 
However, the rm.s. scatter of both the original and the renor- 
malized data set of May et al. is even sUghtly larger than that 
of O'Brian et al., and differences such as in A5395.250 and 
A5487.160 of Mult 1143 with logeFoi = 8.01 and 7.71, re- 
spectively, are also found in their sample. 

Interestingly enough some of the more recent measure- 
ments of the Oxford and Hannover groups seem to produce 
substantially smaller scatter. Whereas cr{\oge) ~ 0.15 for the 
O'Brian et al. and May et al. samples, cr(loge) ^ 0.05 . . . 0.07 
for the Oxford and Hannover lines. Fig. |] shows a marginal dif- 
ference between the two groups, but that depends on a partic- 
ular choice of our models with (A log eFci)oxf-Han = 0.067 
for the HM LTE model and 0.026 for the TH LTE model. Let 
us mention here that line-by-line comparison of /-values of the 
two groups leads to a difference of (Alogg/pcOoxf-Han = 
-0.029 ±0.009. 

In order to evaluate the solar iron abundance we thus de- 
cided to disregard all but the Oxford and Hannover /-values. 




Fig. 8. Logarithmic solar abundances as a function of equiva- 
lent width in mA determined with the HM solar model in LTE 
and ^ = 1.0 kms^^. Top: Oscillator strengths from May et al. 
(1974 filled circles) and from O'Brian et al. ( 1991 , open cir- 
cles). Bottom: /-values from the Hannover group (sources p,q 
of Table 2, filled circles) and from the Oxford group (sources 
e,f,g,h,n,o of Table 2, open circles). The range of ilcr rm.s. 
scatter is indicated by the shading 



Table 3. Solar Fe l abundances based exclusively on the /- 
values of the Oxford and Hannover groups, calculated for dif- 
ferent models of line formation. Note that A log Cg refers to 
Anstee & O'Mara's damping constants. It was chosen so that 
the mean abundances did not depend on equivalent width (see 
left panels in Fig. H). See text for further discussion 





Model 


5 [kms^^] 


AlogCe 


log£Fol,0 





THLTE 


0.85 


-0.12 


7.508 ± 0.080 


1 


NLTE 0+ 


0.85 


-0.23 


7.605 ± 0.087 


2 


NLTE 5+ 


0.85 


-0.10 


7.521 ± 0.089 


3 


NLTE 5- 


0.85 


-0.15 


7.629 ± 0.094 


4 


HMLTE 


1.00 


0.09 


7.574 ± 0.074 


5 


THLTE 


1.00 


-0.14 


7.477 ± 0.070 


6 


NLTE 5+ 


1.00 


-0.12 


7.488 ± 0.075 


7 


NLTE 1 + 


1.00 


-0.13 


7.503 ± 0.077 


8 


NLTE 1 + 


1.00 


-0.16 


7.499 ± 0.075 


9 


NLTE 1/2+ 


1.00 


-0.17 


7.509 ± 0.077 



Unfortunately, this choice reduced our line sample from 391 to 
97 lines. Fig. |^ demonstrates that all of the weak lines in this 
combined sample are from Hannover sources whereas most 
of the strong lines were measured in Oxford. This correlates 
nicely with excitation energies, such that all low-excitation 
lines come from Oxford sources and all high-excitation lines 
are due to Hannover measurements. 

3.2.2. The solar iron abundance 

Irrespective of the choice of the /-values the solar Fe I abun- 
dances as calculated from fitting the solar flux spectrum depend 
sensitively on the model assumptions. Blackwell et al. ( 1995a| ) 
and Grevesse & Sauval (1999) both have reported that the HM 
empirical solar model leads to Fe I abundances systematically 
higher than those obtained from theoretical models or other 
empirical models with a lower temperature in their upper lay- 
ers. This is to be expected under the assumption of LTE since 
the source function then is always Planckian, and the emerging 
intensities in theoretical models will to first order follow the 
temperature stratification. It is, however, not evident for NLTE 
line formation, since there both the source function and the op- 
tical depth scale may deviate from their thermal behaviour 

In Paper I the level populations had been discussed for a 
number of LTE and NLTE population models. It was argued 
there that in most of the NLTE models - at least those with 
non-zero hydrogen collisions - the line source functions were 
very close to thermal, and the differences of line profiles with 
respect to LTE occurred essentially due to parametrization of 
(a) hydrogen collisions and (b) a cutoff energy above which all 
levels were thermalized with respect to the Fe II ground state. 
The latter operation had to be included to simulate the missing 
ionization/recombination channels. The different populations 
are shown in Fig. 6 of Paper I, and as yet we have not been able 
to choose a best case model on the basis of comparison with 
the strong lines only. 

Fig. ^ therefore gives an impression of how the solar Fe I 
abundances obtained from line profile fits based on different 
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Fig. 9. Solar Fe i abundances as a function of W\ (left) and lower-level excitation energy Eiov, (right), calculated for different 
models of line formation. Oscillator strengths are from the Oxford (open circles) and Hannover groups (filled circles). See text 
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Fig. 10. Solar flux spectrum of Fe I 41, A4404.750 A, together 
with three nearly identical synthetic profile fits using the HM 
LTE, the TH LTE and the TH NLTE 1+ models. See text for a 
discussion of strong lines 

LTE and NLTE models with different line-broadening param- 
eters depend on the model assumptions. As mentioned above, 
only the Oxford and Hannover group /-values have been con- 
sidered. With respect to Table 0the models in Fig. ^ are mod- 
ified using the original models 7 and 8 of Table |l] to interpo- 
late corrections of the damping constant so that the resulting 
mean abundances are independent of line strength. As docu- 
mented in Table H these additional corrections are always small. 
Comparing models 7 and 8 in Table || it is evident that the two 
interpolated results do not differ significantly. 

Four characteristic features are displayed in Fig. ^ 

1 . There is still a systematic difference between abundances 
(oscillator strengths) of the Oxford and Hannover groups 
which is seen best in the domain of the turbulence lines 
around 80 to 100 mA. It is also found as a difference 
between lines of low and high excitation. This would be 
even easier to detect if the adjustment of the damping con- 
stants were applied to the individual sets of lines calcu- 
lated from a common base of /-values. For our model 
(4) in Fig. ^ the Oxford data alone then would require a 
damping correction of A log Cq — —0.28, and they would 
lead to the value of logeFci.o — 7.693 ± 0.052, reason- 



ably close to that of Blackwell et al. ( 1995a ). Vice versa, 
LTE in the HM model (4) applied only to the Hannover 
/-values would require AlogCe = +0.13, and result in 
logepoi,© — 7.543 ± 0.070. Our compromise to fit the 
combined data set thus does not at all resolve the long- 
standing discrepancy. It is important to recognize that this 
problem does not seem to depend on the particular LTE or 
NLTE model chosen. The difference between the Oxford 
and Hannover line abundances is only slightly smaller (0. 10 
dex) for the TH LTE model. It is removed here only by ad- 
justment of the damping constants for the individual mod- 
els, the shortcomings of which are hidden in a slightly in- 
creased scatter. 

2. The turbulence lines deviate from both weak and strong 
lines in models (0) to (3), and perhaps in the HM LTE 
model (4) because the mean micro turbulence is relatively 
low. This choice was made in Paper I mostly to model 
a number of the stronger turbulence lines with equivalent 



widths around 100 mA. After having examined a series 
of tests with different values we concluded that a value of 
^ = 1.0 kms^^ produced profile fits of approximately the 
same quality. As is evident from comparing models (0) and 

(5) inFig.|the higher value of ^ tends to improve the uni- 
formity of the abundances. A similar increase would also 
improve the results of the HM model. 

3. Even after having adjusted the strong lines to fit to a com- 
mon mean abundance with the weaker lines it is surpris- 
ing how they lead to systematically lower abundances than 
the sample mean. Part of this difference may be attributed 
to a relatively bad fit of the line wings. Fig. [l^ shows the 
discrepancy between inner and outer wing synthesis. Both 
parts of the profile are of photospheric origin. We note that 
a slightly better fit of the outer wings can be achieved with 
an increase of the iron abundance by ~ 0.03 dex which, 
however, would not remove the trend. Moreover, it would 
destroy the fit of the inner wing to an unacceptable degree. 

4. The run of abundances with excitation energies displays a 
decrease with i^iow for most of the models. As was empha- 
sized by Blackwell et al. (1995a) and Grevesse & Sauval 
(1999) this tendency is relaxed or even removed by intro- 
ducing atmospheric models with lower temperatures in the 
upper photosphere. This trend is confirmed when compar- 
ing the HM and TH LTE models in our analysis. However, 
care must be taken not to confuse it with a similar one pro- 
duced by the dependence upon microturbulence. The cur- 
rent sample of Fe I lines includes quite a number of low- 
excitation lines in the turbulence regime (Mults 1, 2, 3 and 
13), which dominate the least-squares approach in Fig. ^ 
This becomes particularly evident by comparison of the 
LTE model (0) and (5), and by the NLTE models (2) and 

(6) , where the increase of ^ from 0.85 to 1.00 kms^^ re- 
moved most of the energy dependence. 

The solar iron abundance determined by even the most 
careful spectral analysis thus depends on the proper choice of 
both the atmospheric model and the oscillator strengths. While 
Grevesse & Sauval ( 1999| ) claim to have solved the discrepan- 
cies of the long-standing debate on the solar iron abundance 
by introducing their special semi-empirical adjustment to the 
HM atmospheric model, it is only fair to notice that even their 
final data produce an abundance difference with mean values 
of (log£Foi^0)Han = 7.476 ± 0.053, and (logeFei,0)oxf = 
7.514 ± 0.036. What makes this result less useful is the neglect 
of all strong lines. As was shown above it is the strong lines in 
the Oxford sample that - having been adjusted to the weaker 
lines by a corresponding decrease of the damping constants - 
confirm the high solar Fe I abundance claimed by Blackwell 
et al. (1995a). Different from the Kiel-Hannover group the 



Oxford group does not cover the full range of line strengths and 
excitation energies encountered in the solar spectrum. In par- 
ticular the weak lines are missing, for which an analysis would 
allow a direct comparison of the /-value sources without refer- 
ence to the uncertainties of line broadening processes. 

There is no use ignoring the fact that either the oscilla- 
tor strengths currently available are discrepant at a level that 
cannot be explained by laboratory measurement errors alone. 
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or that the solar spectral line identifications are erroneous at 
an equally unacceptable level, or that atmospheric inhomo- 
geneities are much more important for individual lines than 
expected. Let us discuss all three possibilities. 

Much of the different absolute scales of /-values is due to 
the necessary normalization which can be improved; however, 
an individual scatter of lines in a common multiplet is obtained 
even for experimental methods thought to be very accurate. 
As an example let us consider the abundance scatter of lines 
in Mult 1 14. All lines have been measured by the Hannover 
group, and the abundances spread from 7.41 at A5141.739 to 
7.65 at A5049.819 to a value as high as 7.78 for A4924.769 
if the HM LTE model is applied. These are not faint lines for 
which high measurement errors could be accepted; the experi- 
mental error estimates range from 0.04 to 0.07 dex for these 
lines, which transforms to the fact that our abundances lead 
to results that are discrepant on much more than a level. Of 
course, the results may tell us that the hollow-cathode measure- 
ments of A5141.739 are not of the same quality as the other two 
lines which were measured by laser-induced fluorescence, but 
that would invalidate the experimental error estimates. 

Comparison of such multiplet abundance scatter based on 
common source /-values with that already discussed above in- 
dicates that this does not depend very much on the experimental 
methods either, although there may exist still a number of prob- 
lems that are connected with the control of experimental envi- 
ronment parameters as discussed by Holweger et al. (1995). 
Thus we conclude that agreement of mean abundance values 
between different sources of oscillator strengths (often claimed 
for the O' Brian et al. data) is not a significant measure of me- 
thodical accuracies. Taken at face value the r.m.s. scatter of 
abundances obtained from a single set of oscillator strengths 
such as that of O' Brian et al. is a measure of the accuracy of 
the mean solar Fe I abundance that can be reached with these 
data. In fact the accuracy is then even less due to blends and 
other problems referring to the profile fits, and to the ambigui- 
ties of atmospheric modelling. 

There exists a number of lines in the iron spectrum that 
could be misidentified in that the spectral features could be 
blends that are not only unresolved but also fall within a few 
mA of the same center wavelength. As with other undetected 
blends such profiles will be fitted with too large abundances. 
This should produce abundance distributions that are system- 
atically shifted to the high-abundance side, something that is 
not detected in the results. To reduce the dominating intrinsic 
abundance scatter to reasonable amounts it would mean that 
more than half of the lines would have to be corrected for such 
blend or identification problems, a situation that seems highly 
unlikely. We note that many blend problems of the kind pro- 
ducing too large fit abundances are avoided by our profile fit- 
ting method which allows an exchange of certain fit parame- 
ters such as abundance, microturbulence or damping parame- 
ters only within a narrow region. In such cases the profile fit 
procedure always tends to produce higher abundances. 

Our discussion of line broadening in subsection 3.1.2 and 
Fig. H has shown that the true abundance differences resulting 
from line formation in plane -parallel and in hydrodynamic at- 
mospheres are quite small. They are even negligible taking into 



account the large abundance differences that appear between 
sets of different /-values. The mere change of atmospheric 
models affects the mean abundance but not the r.m.s. scatter 
as can be found in Table ||, and it is obvious that changing the 
microturbulence has a greater influence on such results. Thus it 
is doubtful if any other atmospheric model could significantly 
reduce the abundance scatter. 

Our results then indicate that it is the atomic data, in par- 
ticular the oscillator strengths, that presently do not allow the 
determination of the solar Fel abundance with an accuracy 
better than ^ 0.1 dex. Based on the most reliable sets of /- 
values (Oxford and Hannover data) and on the model produc- 
ing the smallest overall dependence on excitation energy (TH 
NLTE 1/2H-) we find a value of logeFei,© = 7.509 ± 0.075 
with no dependence on line strength but a small residual gra- 
dient with energy, Aloge/Ai^iow.cV = —0.005. In view of 
the differences between the Oxford and Hannover /-values it 
is important to notice that this value is only 0.02 dex above 
that obtained from the Hannover data alone, while it is 0.09 
dex below the pure Oxford value. This apparent contradiction 
is resolved by inspection of the corresponding energy depen- 
dence of the respective sources. Whereas the Hannover re- 
sults show no energy gradient, the Oxford data - after hav- 
ing adjusted the damping constants to remove a line strength 
trend - keep a strong gradient with excitation energy for which 
Aloge/AiJiow.cV = 0.034. The last three models in Fig. ^ 
show only a small residual energy dependence of the Fe I abun- 
dances ranging from A log e / AiJiow.cV = —0.0094 for the TH 
NLTE5H- model to Aloge/A£;iow,oV = -0.0054 for the TH 
NLTE1/2H- model. 

The above results are to be understood as a clear report 
of our failure to solve the photospheric solar Fe I abundance 
problem if more than the Hannover data set were involved. 
Using this data set alone with the HM LTE model, a micro- 
turbulence of 1.05 kms^^ together with damping corrections 
AlogCe — 0.11 (above the Anstee & O'Mara damping con- 
stants) yields log epci,© — 7.535 ± 0.070. The energy gradient 
for that result is Aloge/Ai?iow,GV = —0.008. The overall best 
NLTE model (TH NLTE 1/2h-) appHed to the Hannover data 
alone leads to logepci,© = 7.480 ± 0.072 with no dependence 
on energy. 

4. Conclusions 

The choice of a particular model to determine the solar Pel 
line formation with a valid parametrization of the atomic colli- 
sions is not possible even when including the weak solar lines. 
Arguments referring only to the solar abundance problem with 
or without inclusion of the Fe II lines are not conclusive since 
both sets of /-values (Fe I and Fe ll) are far from producing ho- 
mogeneous results. One marginal result is that the models of 
Paper I with their low microturbulence are no longer compet- 
itive because they all display a relatively strong gradient with 
excitation energy (see Fig. ^). This does no longer appear when 
increasing the microturbulence from ^ — 0.85 kms^^ to 1.00 
km s^^ as in our present models 5 to 9. All the TH models are 
roughly compatible with meteoritic abundance. Small correc- 
tions for dynamic line formation such as suggested by com- 
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parison with hydrodynamic results of Asplund et al. (2000) in 
section 3.1.2 are of the order of —0.03, which would bring the 
solar abundance to a value slightly below that of the carbona- 
ceous chondrites. 

The quality of individual line fits are significantly differ- 
ent for the HM and TH model atmospheres only for the cores 
of strong lines. In Paper I this was demonstrated for a num- 
ber of lines of various excitation energies. The line center flux 
reflects essentially the different temperatures in the upper pho- 
tosphere with a 150 .. . 200 K difference predicting AF ^ 4% 
as observed. However, these differences vanish when a com- 
promise is accepted for a profile fit of the inner wings (see Fig. 

allowing the synthetic profile to fall below the observed 
flux by a small amount. The evaluation of profile fits thus has 
changed marginally as compared with Paper I. For the weaker 
lines Figs. Q and g document the independence of fit quality 
from the model atmosphere if abundances and macroturbulence 
velocities are adjusted accordingly. 

The selection of a particular atmospheric/atomic model on 
the grounds of profile synthesis of the solar Fe I flux spectrum 
is therefore still somewhat ambiguous. This would be differ- 
ent if the abundance determinations were of higher quality. For 
differential analyses of stellar spectra it is obvious that our at- 
mospheric model can be only one of the TH models because 
only they allow a physically consistent change of parameters 
such as Teff, log g or [Fe/H]. Since strong lines in the solar 
spectrum reduce to weak or turbulence lines in stars of low 
metal abundance, it is most important to install a unique recipe 
for the determination of the damping parameter. This can be 
done with reference to Tabled where a good mean value for the 
correction would be AlogCg = —0.15. We should, however, 
bear in mind that this deviation from the Anstee & O'Mara 
results is essentially necessary to coiTect the strong fines with 
/-values from the Oxford group. The eiTor introduced to differ- 
ential abundance determinations in metal-poor stars thus will 
have to include a systematic uncertainty of ^ 0.04 dex due to 
inconsistencies in the interpretation of the solar lines. 

CuiTent investigations of a small number of reference stars 
with different iron abundances will have to show how to se- 
lect a common NLTE model that fits the Fe Il/Fe I ionization 
equilibria of all stars. 
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